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ABSTRACT 


Reflection,  transmission  and  absorption  of  solar  radiation  by  -fotir>  y ) 

, o.  of'"*  y 

cloud  types  (low  cloudy  middle  JcT©ujL/,h1gh  cloudiland  stratus)  are  com- 
puted as  functions  of  the  solar  zenith  angle  and  cloud  liquid  water/ice 

-r 

content^  ^The  reflection,  transmission  and  emission  of  infrared  radia- 
tion by  cirrus  clouds  are  calculated  as  functions  of  th£  cloud  ice  con- 
tent. The  plane-parallel  radiative  transfer  program  employed  Is  based 
on  the  discrete-ordinate  method  with  applications  to  Inhomogeneous 
atmospheres  covering  the  entire  solar  and  Infrared  spectrums  and  taking 
into  account  the  gaseous  absorption  In  scattering  atmospheres. 

— A.  ^ 

The  resulting  values  of  the  solar  radiative  properties  of  clouds 
are  fitted  with  known  mathematical  functions  involving  the  solar  zenith 
angle  and  cloud  liquid  water/ice  content  as  variables.  Parameterized 
equations  for  the  infrared  flux  reflectivity,  transmissivity,  and 
emissivity  of  cirrus  clouds  are  *Tso  presented  as  functions  of  the 
cloud  ice  content.  Effects  of  the  atmospheric  profile  are  discussed 
and  the^effects  of  surface  reflectivity  on  the  solar  radiative  proper- 
ties  of  clouds  are  parameterized  In  terms  of  the  water  vapor  absorp- 
tivity below  the  cloud,  ground  reflection  and  average  cloud  reflection. 

The  parameterized  empirical -theoretical  functions  for  the  solar 

r < file* 

and  Infrared  radiative  properties  of  clouds  are ^t+Hzedyn  a global 
radiation  budget  model.  The  resulting  radiation  balance  of  the  earthj_^/ 
atmosphere  system  matched  ’Ste^satell 1 te  observations  of  planetary  — 
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albedo,  absorption  of  solar  radiation,  and  upwelling  Infrared  radiation 
more  closely  than  did  earlier  computed  models.  The  planetary  albedo 
was  calculated  to  be  0.31. 
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CHAPTER  1 


INTRODUCTION 

The  most  important  regulators  of  the  radiation  balance  are  clouds, 
which  regularly  occupy  at  least  fifty  percent  of  the  sky  on  a global 
scale.  Clouds  absorb  and  scatter  the  Incoming  solar  radiation  and  emit 
thermal  Infrared  radiation.  The  amount  of  energy  absorbed  and/or  emit- 
ted represents  one  of  the  prime  sources  determining  the  stability  of 
cloud  layers  and  Is  further  associated  with  the  general  circulation  of 
the  atmosphere.  Moreover,  the  albedo  of  the  earth-atmosphere  system 
depends  crucially  on  such  parameters  as  cloud  cover,  cloud  composition 
and  cloud  structure. 

Reflection,  transmission  and  absorption  of  solar  radiation  by  var- 
ious cloud  types  imbedded  In  atmospheres  containing  water  vapor  was 
extensively  investigated  by  Llou  (1976),  and  recently,  Llou  et  al. 
(1978)  extended  the  radiative  transfer  program  to  include  the  absorp- 
tion contribution  of  oxygen,  ozone  and  carbon  dioxide  and  the  absorp- 
tion and  scattering  contributions  by  aerosols  typical  of  those  In  clear 
atmospheres.  The  radiation  program,  as  utilized  In  this  research,  has 
taken  into  account  the  inhomogeneity  of  the  plane-parallel  atmosphere, 
the  wavelength  dependence  of  solar  radiation  and  the  gaseous  absorption 
In  scattering  layers.  Calculations  have  been  carried  out  to  cover  the 
entire  solar  spectrum  on  a band-by-band  basis. 

Although  the  radiation  characteristics  of  clouds,  whose  horizontal 


2 


extents  are  much  greater  than  their  vertical  depths,  can  be  evaluated 
to  yield  good  accuracies  by  means  of  current  existing  transfer  methods, 
the  computational  time  requirement  and  the  effort  involvement  make  the 
comprehensive  transfer  program  impracticable  for  studies  of  the  weather, 
climate  and  radiation  balance  of  the  earth-atmosphere  system.  Owing  to 
the  intricacy  of  the  cloud  interaction  with  the  solar  radiation  field 
of  the  atmosphere,  a set  of  prescribed  values  for  reflection  and  trans- 
mission values  has  normally  been  used  In  the  study  of  the  thermodynamic 
properties  of  the  atmosphere.  Manabe  (1975),  for  example,  assumed  that 
the  reflection  of  high,  middle  and  low  clouds  were  20*,  48%  and  89%, 
respectively,  regardless  of  the  cloud  composition  and  solar  zenith 
angle. 

It  Is  the  purpose  of  this  research  to  provide  empirical-theoretical 
equations  for  the  reflection,  transmission  and  absorption  of  solar  radi- 
ation by  clouds  of  different  type  as  functions  of  the  solar  zenith 
angle  and  cloud  liquid  water/ice  content.  The  mathematical  known  func- 
tions are  derived  by  means  of  numerical  fits  of  the  precalculated  re- 
flection, transmission  and  absorption  values  for  a black  surface. 
Parameterization  of  the  effects  of  surface  reflection  and  absorption 
by  water  vapor  between  the  cloud  and  the  earth's  surface  is  given  In 
terms  of  the  water  vapor  absorptivity,  ground  reflection  and  average 
cloud  reflection.  Parameterized  equations  for  the  Infrared  ilux  reflec- 
tivity, transmissivity  and  emlsslvlty  of  a cirrus  cloud  are  also  pre- 
sented as  functions  of  the  cloud  vertical  Ice  content.  In  addition, 
these  parameterized  equations  are  used  In  conjunction  with  the  radiative 
transfer  model  program  to  determine  the  radiation  balance  of  the  earth- 
atmosphere  system. 
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CHAPTER  2 


RADIATIVE  TRANSFER  IN  A CLOUDY  ATMOSPHERE 

The  basic  radiation  scheme  used  is  based  on  the  plane-parallel 
model  recently  developed  by  Liou  et  al.  (1978),  Freeman  and  Liou  (1978), 
and  Roewe  and  Liou  (1978).  The  model  utilizes  the  discrete-ordinate 
method  for  radiative  transfer,  originally  introduced  by  Chandrasekhar 
(1950).  The  method  has  been  found  quite  satisfactory  because  it  allows 
for  the  solutions  of  the  integral -differential  transfer  equation  to  be 
explicitly  derived  and  can  easily  take  into  account  the  distribution, 
thickness  and  types  of  clouds  and  aerosols,  making  the  scheme  well 
suited  to  this  study. 


2.1  Method  of  Radiative  Transfer 

The  transfer  equations  describing  the  radiation  field  for  the 
azimuthal  independant  diffuse  intensity  I may  be  written 

u » I(t,u)  - y-  7 P(y,y‘)I(y,y' )dy' 

T F0p(w*wo)exp(~T/y0)  (solar) 

(1  -ai0)Bv(T)  (Infrared)  (1) 

where  r represents  the  optical  depth,  y and  yQ  the  cosine  of  the  emer- 
gent and  solar  zenith  angles,  respectively,  <Lo  the  single  scattering 
albedo,  FQ  the  incident  solar  flux,  T the  temperature;,  P the  normalized 


phase  function  and  Bv  the  Planck  function.  All  parameters,  except  T, 
y and  yQ  are  wavelength  or  wavenumber  dependent. 

Employing  the  discrete-ordinate  method  for  radiative  transfer, 
the  solutions  of  the  transfer  equations  as  shown  by  Liou  (1973)  are 
given  by 


-k,x 


Kt.Vj)  = J L (^1  )e  ^ 

Z(yi)expC-T/u0) 


+ < 


(solar) 

(Infrared)  (2) 

where  i and  j (-n,n)  denote  the  discrete  streams  (positive  upward  and 


BV<T> 


negative  downward)  used  in  approximating  the  basic  transfer  equation, 
<j)j  and  kj  represent  the  eigenfunction  and  eigenvalue,  respectively, 
the  Z-function  is  associated  with  Chandrasekhar's  H-functlon,  the 
single  scattering  albedo  and  the  expanded  phase  function  In  terms  of 
Legendre  polynomials,  and  L/s  are  coefficients  determined  from  the 
boundary  conditions. 

The  solutions  of  the  radiative  transfer  equations  given  by  Eq. (2) 
are  applicable  to  a homogeneous  and  isothermal  layer.  In  order  to 
apply  such  solutions  to  Inhomogeneous  atmospheres,  the  atmosphere  is 
divided  into  a number  of  sub-layers  each  of  which  may  be  considered 
homogeneous  and  Isothermal.  At  the  top  of  the  atmosphere  there  is  no 
downward  diffuse  tntensity  so  that 

II (o ,-u1 ) » 0.  (3) 


Between  the  sub-layers,  the  Intensities  from  all  directions  must  be 


continuous.  Thus, 


= It+l^Tru1^’  1 3 ]»2 N'1  (4) 

where  N Is  the  total  number  of  sub-layers  and  represents  the  opti- 
cal depth  from  the  top  of  the  atmosphere  to  the  layer  l.  At  the  bot- 
tom of  the  atmosphere,  the  upward  intensity  is  given  by 


VvV  ' 


A r n 

— |2*.Z  *i'“iI(V*u1)  + n.0F0  exp(-rN/v0) 

(so 


ar) 


B (T  ) 
v s 


(Infrared) 


(5) 


where  A$  represents  the  surface  albedo,  ai  and  are  Gauss  weighting 
factors  and  zeros  of  the  Legendre  polynomials,  respectively,  n the 
number  of  discrete  streams  and  Ts  surface  temperature. 

Upon  Inserting  the  intensity  solution  into  Eqs.  (3)-(5),  a set  of 
linear  equations  Is  obtained  from  which  the  unknown  Lj  for  each  sub- 
layer may  be  determined  by  means  of  a matrix  inversion  technique. 

■After  the  constants  of  proportionality  have  been  derived,  the  intensity 
distribution  at  any  level  in  the  atmosphere  may  be  evaluated. 

The  upward  and  downward  fluxes  are  then  given  by 


F+(t)  ■ 2ir  ■ J a^u^Hx.u^) 

F4'(t)  - -2m  l a.p.I(x,-p.)  - S(x,-p  ) 
1-1  1 1 1 0 


where 


(6) 


' ' 

■MM 


*■  ~ 


6 


T 


hni0F0  exp(-T/yQ)  (Solar) 

s(T>-u0)  = \ 

i 0 (Infrared) 

Finally,  the  net  flux  for  a given  layer  may  be  computed  from 
F(t)  = F+(t)  + F+(t)  . 


(7) 


(8) 


n 

ii 


2.2  Solar  Spectrum 

The  primary  gaseous  absorbers  considered  in  the  solar  spectrum  are 
water  vapor,  ozone,  carbon  dioxide  and  molecular  oxygen.  In  addition, 
scattering  and  absorption  by  clouds  and  aerosols  were  taken  Into  ac- 
count. The  range  of  the  solar  spectrum  considered  in  this  study  was 
that  from  0.2  to  3.4  micrometers  as  plotted  In  Figure  1,  where  the 
inner  curve  Is  the  solar  spectrum  at  the  bottom  of  the  atmosphere,  and 
the  shaded  areas  are  absorption  bands.  The  particular  bands  chosen, 
the  primary  absorber  In  each  band  and  the  fractional  solar  flux  In  each 
band  (from  Thekaekara,  1974)  are  given  In  Table  1. 

Table  1.  Solar  Bands,  Absorbers  and  Fractional  Solar  Flux 

BAND  (urn)  0.3  0.5  0.7  0.94  1.10  1.38  1.87  2.70  3.20 

ABSORBER  03  03  02  H20  H20  HgO  C02/H20  HgO 

FRACTIONAL 

SOLAR  FLUX  .0872  .2696  .2034  .1347  .0892  .1021  .0622  .0300  .0216 


2.3  Infrared  Spectrum 

The  primary  gaseous  absorbers  considered  by  this  research  In  the 
infrared  spectrum  are  water  vapor,  carbon  dioxide  and  ozone.  Figure  2 
displays  the  infrared  spectrum  as  obtained  by  satellite,  superimposed 


TFi 

r'r*  • 


WAVENUMBER  (cm 


9 

on  plots  of  the  spectral  distribution  of  thermal  emission  for  several 
earth-like  temperatures.  The  infrared  absorption  bands  used  in  this 
study  are  labeled  on  the  plot  and  listed  in  Table  2. 


Table  2.  Infrared  Absorption  Bands  and  Fractional  Infrared  Flux 


Absorption  Band 

Wavenumber  Range 

Fractional  Infrared 

Flux 

1. 

^0  rotational 

40-900 

0.615 

2. 

15  pm  COg 

582-752 

0.150 

3. 

HgO  continuum 

800-1200 

0.160 

4. 

9.6  ym  0^ 

1000-1065 

0.024 

5. 

6.3  pm  HgO 

1200-2200 

0.051 

2.4  The  Model  Atmospheres 

Five  model  atmospheres  (McClatchey  et  al.,  1971)  were  used  in  this 
study.  The  atmospheres  and  the  latitude  belts  to  which  they  apply  are 
listed  in  Table  3.  Each  atmosphere  includes  vertical  profiles  of  pres- 
sure, temperature,  density,  water  vapor  and  ozone.  Figures  3 and  4 
present  the  atmospheric  profiles  of  temperature  and  ozone  density. 
Figure  5 displays  the  water  vapor  distributions  and  Figure  6 the 


Table  3.  Five  Model  Atmospheres 


Atmosphere 

Latitude  Belts 

Subarctic  Winter  and  Summer 

90°N-60°N  and  90°S-60°S 

Midlatitude  Winter  and  Summer 

60°N-30°N  and  60°S-30°S 

Tropical 

30°N-30°S 

TROPICAL 


WATER  VAPOR  DENSITY  (gm  m s)  (cm's) 

Figure  5.  Mater  vapor  density  profiles  in  the  Figure  6.  Aerosol  and  molecular 

model  atmospheres.  distributions. 


] 


12 

profiles  of  aerosol  and  molecular  densities.  Table  4 lists  the  con- 
centrations of  the  uniformly  mixed  gases,  carbon  dioxide  and  molecular 
oxygen. 


Table  4.  Concentrations  of  the  Uniformly  Mixed  Gases 


Constituent 

Molecular  Wt. 

ppm  by  Volume 

_2 

g cm  /mb 

Air 

28.97 

106 

1.02 

co2 

44 

330 

5.11  X 10'4 

CVJ 

O 

32 

2.095  X 105 

0.236 

2.5  Physical  Properties  of  the  Cloud  Types 

The  effects  of  clouds  on  the  transfer  of  solar  and  infrared  radia- 
tion are  determined  by  the  particle  phase  (liquid  water  or  ice),  con- 
centration, size  and  size  distribution.  These  factors  combine  to  de- 
termine the  single  scattering  albedo  and  phase  function. 

Four  different  cloud  types  were  utilized  in  this  research.  The 
cloud  height  values  (London,  1957)  and  other  cloud  parameters  used  are 
listed  in  Table  5.  The  cirrus  clouds  were  considered  to  be  composed 


Table  5.  Cloud  Types  and  Physical  Properties 


Cloud  Type 

Concentration 

(an-3) 

Drop  Radius 
(pm) 

Water/ Ice  Content 
(g  m'3) 

Cloud  Base 
(km) 

Cumulus 

300 

0-40 

0.33 

1.8 

A1 tostratus 

450 

0-13 

0.24 

4.1 

Stratus 

178 

0-40 

0.78 

1.4 

Cirrus 

0.1 

0.05 

10.2 

I 

I 

I 

I 

I 

I 

I 

r 

i 
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i 
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exclusively  of  ice  cylinders,  randomly  oriented  in  the  horizontal  plane, 
with  a mean  length  of  200  pm  and  a mean  radius  of  30  pm.  The  drop-size 
distributions  for  the  cumulus,  altostratus  and  stratus  clouds  were 
based  on  the  observations  of  Battan  and  Reitan  (1957),  Diem  (1948)  and 
Singleton  and  Smith  (1960),  respectively.  These  size  distributions  are 
shown  in  Figure  7,  the  curves  being  presented  in  such  a way  that  the 
integral  of  the  area  under  each  curve  Is  equal  to  the  total  particle 
concentration  indicated  in  Table  5.  All  clouds,  for  the  purposes  of 
the  radiative  transfer  calculations,  are  considered  to  be  of  infinite 
horizontal  extent  and  plane  parallel.  Within  the  infrared  spectrum  all 
clouds  except  cirrus  clouds  were  considered  to  be  black  bodies. 

2.6  Radiative  Properties  of  Clouds 

As  described  previously,  the  radiation  model  calculates  the  trans- 
fer of  solar  radiation  in  a series  of  absorption  bands.  To  obtain  the 
reflection,  transmission  and  absorption  of  the  cloud  for  the  entire 
solar  spectrum,  proper  summation  over  the  fluxes  in  each  sub-spectral 
region  weighted  by  the  appropriate  percentage  of  solar  flux  is  required. 

The  reflection  r may  then  be  defined  as  the  ratio  of  the  reflected 
flux  to  the  incident  solar  flux  normal  to  the  cloud  top.  Hence,  the 
reflection  of  a cloud  l?yer  for  the  entire  solar  spectrum  is  given  by 

r-$foX>V<Tt>][vsJ-  <9> 

where  f^  denotes  the  amount  of  solar  flux  In  the  ith  spectral  band,  SQ 
the  solar  constant  and  rt  the  optical  depth  at  the  cloud  top. 

Similar  to  the  above  definition,  the  transmission  t of  a cloud 


' 
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layer  may  be  defined  as  the  ratio  of  the  downward  flux  at  the  cloud 
base  to  the  incident  solar  flux  normal  to  the  cloud  top.  Thus,  the 
transmission  of  a cloud  layer  for  the  entire  solar  spectrum  is  given  by 

t=^*K>/F>t>][VSo]  • HO) 

where  denotes  the  optical  depth  at  the  cloud  base. 

Moreover,  the  absorption  a of  solar  flux  within  a cloud  layer  can 
be  evaluated  from  the  net  flux  divergence  between  the  top  and  bottom 
of  the  cloud.  Therefore,  the  total  absorption  within  a cloud  layer  for 
the  entire  solar  spectrum  is  given  by 


The  radiation  model  also  calculates  the  transfer  of  infrared  radi- 
ation in  a series  of  absorption  bands.  Reflectivity,  transmissivity 
and  emissivity  of  the  cloud  for  the  entire  infrared  spectrum  may  also 
be  defined  in  terms  of  summations  over  the  fluxes  in  each  sub-spectral 
region. 

The  flux  reflectivity  rp  of  the  cirrus  cloud  is  defined  as  the 
ratio  of  the  downward  flux  at  the  cloud  base  to  the  incident  infrared 
flux  at  the  cloud  base.  Therefore,  the  flux  reflectivity  of  the  cirrus 
cloud  for  the  entire  Infrared  spectrum  Is  given  by 

('2> 

f h 

where  f^  Is  the  amount  of  infrared  flux  in  the  i spectral  band  and 
the  total  Infrared  flux. 

• o 
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The  flux  transmissivity  tp  is  defined  as  the  ratio  of  the  upward 
flux  at  the  cloud  top  to  the  incident  infrared  flux  at  the  cloud  base. 
Thus,  the  flux  transmissivity  of  the  cirrus  cloud  for  the  entire  infra- 
red spectrum  is  given  by 

tF*I  (13) 

Furthermore,  the  flux  emissivity  of  the  cloud  layer  is  defined 
in  terms  of  the  emission  from  the  top  and  the  base  of  the  cloud.  The 
flux  emissivity  at  the  cloud  top  £pt  is  defined  as  the  ratio  of  the 
emitted  flux  from  the  cloud  top  to  the  black  body  flux  for  the  cloud 
top  temperature.  Similarly  the  flux  emissivity  at  the  cloud  base  epb 
is  defined  as  the  ratio  of  the  emitted  flux  from  the  cloud  base  to  the 
black  body  flux  for  the  temperature  of  the  cloud  base.  Hence,  the  flux 
emissivities  of  a cirrus  cloud  for  the  entire  Infrared  spectrum  are 


given  by 

[vRo] 

(14a) 

EFb  - | [VRo] 

(14b) 

where  Tt  and  Tfa  are  the  temperatures  at  the  cloud  top  and  base,  respec- 
tively. 

It  is  Important  to  note  that  the  flux  transmissivity  tp  and  the 
flux  reflectivity  rp  represent  quantities  without  Including  the  cloud 
emission.  These  values  are  evaluated  by  letting  B^-(T)  * 0 within  the 
cloud  layer.  On  the  other  hand,  the  flux  emissivities  eFt  and  Cpb  are 
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calculated  using  only  the  emission  due  to  the  cloud  layer.  These 
quantities  are  evaluated  by  letting  F^t^)  and  F+(xt)  equal  zero. 
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CHAPTER  3 

PARAMETERIZATION  OF  REFLECTION,  TRANSMISSION  AND 
ABSORPTION  OF  SOLAR  RADIATION  BY  CLOUD  LAYERS 

The  radiative  transfer  model  discussed  in  the  previous  chapter 
was  used  to  compute  reflection,  transmission  and  absorption  of  solar 
radiation  by  each  of  the  four  cloud  types.  The  following  sections  of 
this  chapter  include  discussions  of  the  effects  of  cloud  thickness, 
atmospheric  profile  and  cloud  type  on  the  solar  radiative  properties 
of  clouds,  and  the  parameterization  of  these  properties. 

3.1  Radiative  Transfer  Model  Computations  in  the  Solar  Spectrum 
Computations  were  made  for  various  cloud  thicknesses,  holding  the 

cloud  base  at  a constant  height  for  each  case.  The  vertical  liquid 
water  content  W of  each  cloud  thickness  is  given  by 

W = w A z,  (15) 

where  w is  the  measured  water/ice  content  given  in  Table  5 and  Az  Is 
the  geometrical  cloud  thickness.  Table  6 lists  the  range  of  cloud 
thicknesses  and  vertical  liquid  water/ice  contents  used  in  the  radia- 
tive transfer  calculations  of  this  study. 

3.2  Effects  of  Atmospheric  Profile,  Cloud  Thickness  and  Cloud  Type 
The  resulting  reflection,  transmission  and  absorption  by  the  cloud 

layer  were  initially  obtained  for  cumulus  and  stratus  clouds  in 
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Table  6.  Range  of  Cloud  Thicknesses  and  Vertical 
Liquid  Water/Ice  Contents 


Cloud  Type 

Thickness 

(km) 

Vertical  Liquid  Water/Ice 
Content 

(g-nr2) 

Cumulus 

0.15  - 2.25 

49.5  - 742.5 

Altos tratus 

0.10  - 1.50 

24.0  - 360.0 

Stratus 

0.50  - 0.75 

39.0  - 585.0 

Cirrus 

0.10  - 2.90 

5.2  - 150.4 

tropical,  midlatitude  winter  and  subarctic  winter  atmospheres  with  a 
zero  surface  albedo.  Figures  8 and  9 show  reflection,  transmission 
and  absorption  of  solar  radiation  by  the  cloud  layer  as  functions  of 

the  cosine  of  the  solar  zenith  angle  for  four  thicknesses  of  a cumulus 

and  stratus  cloud,  respectively.  The  values  of  W-j , and  are 

50,  150,  250  and  450  g-m"^,  respectively. 

The  reflection, transmission  and  absorption  values  for  both  cloud 
types  show  no  significant  dependence  upon  the  atmospheric  profile.  In 

each  case,  the  variance  is  less  than  2%.  The  cloud  thickness  or  the 

vertical  liquid  water  content  is,  however,  a significant  parameter 
whose  change  greatly  affects  the  reflection,  transmission  and  absorp- 
tion of  solar  radiation.  For  a change  in  vertical  liquid  water  content 

mO 

of  400  g-m  (0.5  km  change  In  cloud  thickness),  reflection  varies  by 
as  much  as  43%,  transmission  56%  and  absorption  12%  when  the  sun  is 
overhead  the  stratus  cloud. 

Figure  10  depicts  the  solar  radiative  properties  for  several  cases 
of  cumulus  and  stratus  clouds  having  the  same  vertical  liquid  water 
content.  The  variations  in  reflection,  transmission  and  absorption 


radiative  properties  of  a stratus  cloud  In  tropical  and  midlatitude  winter  atmospheres 
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for  the  two  cloud  types  are  due  to  cloud  compositions,  primarily  the 
particle  size  distribution.  As  shown  in  Figure  7,  the  distributions 
are  such  that  the  cumulus  has  a significantly  greater  number  density 
of  particles  of  less  than  12  pm  diameter,  while  the  stratus  has  a 
slightly  greater  number  density  of  particles  of  greater  than  12  ym 
diameter.  From  Figure  10,  it  is  evident  that  the  cloud  reflection, 
transmission  and  absorption  depend  significantly  upon  the  cloud  type, 
with  some  variations  exceeding  15%  when  the  sun  is  overhead. 

3.3  Parameterization  Procedures 

Based  on  the  results  presented  in  the  previous  section,  calcula- 
tions  were  made  for  each  of  the  four  cloud  types  in  a tropical  atmos- 
phere with  a zero  surface  albedo.  The  cloud  reflection,  transmission 
and  absorption  values  were  computed  as  functions  of  the  cosine  of 
solar  zenith  angle  yQ  and  vertical  liquid  water/ice  content  W.  The 
solar  radiative  properties  were  computed  for  six  values  of  yQ  (0.01, 
0.2,  0.4,  0.6,  0.8,  1.0)  and  fifteen  thicknesses  of  each  cloud  type 
over  the  ranges  listed  in  Table  6.  Figures  11  and  12  are  three- 
dimensional  plots  of  the  resulting  reflection,  transmission  and  absorp- 
tion of  solar  radiation  by  cumulus  and  stratus  clouds,  respectively. 

The  reflection,  transmission  and  absorption  surfaces  for  each  of 
the  cloud  types  were  then  fit  to  an  approximating  polynomial  equation 
using  a bivariate  regression  analysis.  The  coefficients  of  each  of 
the  approximating  predictors  within  the  polynomial  equation  were  cal- 
culated such  that  the  sum  of  the  square  of  the  differences  between  the 
actual  values  of  the  surface  and  values  computed  from  the  polynomial 
equation  was  a minimum. 
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Figure  11.  Solar  radiative  properties  of  a cumulus  cloud 
as  functions  of  vertical  liquid  water  content 
and  cosine  of  solar  zenith  angle. 
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The  following  polynomial  equation  was  found  to  yield  the  best  fit 
of  the  data: 


S<VW)  = l l bijyo 

0 i=0  j=0  13  0 


V 


(16) 


where  S(uQ,W)  denotes  the  cloud  reflection,  transmission  or  absorption, 

O O 

b.j  are  the  predictor  coefficients  and  W is  in  units  of  10  g-m  . 

Tables  7 through  10  list  the  coefficients  which  were  computed  for 
a cumulus,  altostratus,  stratus  and  cirrus  cloud,  respectively.  The 
root  mean  square  error  of  the  approximating  functions  is  less  than  2% 
for  reflection  and  transmission  and  less  than  1%  for  absorption. 
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Table  7.  Approximating  predictor  coefficients,  b.., 

® J 

for  a cumulus  cloud 


REFLECTION 

RMS 

Error  =1.5% 

N 

i o 

1 

2 

3 

0 

.51771E+00 

. 18270E+00 

-.35851E-01 

.23478E-02 

1 

.72387E+00 

- . 271 09E+00 

.51 303E-01 

-.34428E-02 

2 

-. 21499E+01 

.75225E+00 

-.14588E+00 

.97306E-02 

3 

. 12157E+01 

-.40120E+00 

.78145E-01 

-. 52224E-02 

TRANSMISSION  RMS  Error  =1.7% 


1 0 

1 

2 

3 

0 

. 43642E+00 

-.22727E+00 

. 44702E-01 

-.29117E-02 

1 

-.66319E+00 

.35415E+00 

- . 70708E-01 

.46403E-02 

2 

.19014E+00 

-. 10256E+01 

.20574E+00 

- . 1 3534E-01 

3 

10157E+01 

. 54533E+00 

-.  10912E+00 

.71667E-02 

0 

1 

2 

3 


•66103E-01 
. 79484E-01 
. 29293E+00 
. 22077E+00 


ABSORPTION  RMS  Error  = 0.6% 


1 


.34702E-01 
- . 7361 4E-01 
.25047E+00 
-. 13272E+00 


2 


-.69662E-02 
. 17572E-01 
- . 55398E-01 
.28740E-01 


3 


.44253E-03 
-.10788E-02 
. 351 39E-02 
-. 17990E-02 


i 


Table  8.  Approximating  predictor  coefficients,  b^, 

for  an  altostratus  cloud  & 


REFLECTION 

RMS 

Error  *1.3% 

\ 

0 

1 

2 

3 

0 

.64457E+00 

.28785E+00 

-.12069E+00 

. 16673E-01 

1 

•44388E+00 

-.45496E+00 

. 18695E+00 

-.26672E-01 

2 

14875E+01 

. 1 3141E+01 

-.54825E+00 

. 771 21 E-01 

3 

.85585E+00 

-.71346E+00 

. 29823E+00 

- . 41 922E-01 

TRANSMISSION 

RMS  Error  =1.5% 

0 

Y — - 

1 

2 3 

0 . 30735E+00  -.35373E+00  .14797E+00  -.20208E-01 

1 - . 44042E+00  . 501 88E+00  -.20975E+00  .28630E-01 

2 . 13537E+01  -.15738E+01  .66248E+00  -.90739E-01 

3 -. 73396E+0J  .85256E+00  -.35858E+00  .49077 E-01 


ABSORPTION 

RMS 

Error  = 0.4% 

0 

1 

2 

3 

0 

.70594E-01 

.40901 E-01 

- . 1 6982E-01 

.21 345E-02 

1 

-.28640E-01 

- . 1 891 7E-01 

•111 67E-01 

-.37025E-03 

2 

. 1 9500E+00 

. 1 9298E+00 

-.86520E-01 

.98269E-02 

3 

-.15133E+00 

-.10728E+00 

.47143E-01 

- . 53524E-02 
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Table  9.  Approximating  predictor  coefficients,  b^, 
for  a stratus  cloud 


REFLECTION  RMS  Error  >1.5* 


0 

1 

2 

3 

0 

. 38094E+00 

. 26622E+00 

-.62460E-01 

.49964E-02 

1 

. 88882E+00 

-.39997E+00 

. 951 33E-01 

-.76756E-02 

2 

-.26293E+01 

.11593E+01 

- . 28004E+00 

. 22779E-01 

3 

•14876E+01 

-.64913E+00 

.15951 E+00 

- . 1 31 06E-01 

i 

TRANSMISSION 

RMS 

Error  =1.5* 

\ 

0 

1 

2 

3 

0 

.57723E+00 

- . 32665E+00 

.761 60E-01 

-.60763E-02 

1 

- . 88652E+00 

. 52284E+00 

- . 1 2434E+00 

. 10020E-01 

2 

.25788E+01 

-.15559E+01 

.37421 E+00 

- . 30349E-01 

3 

- . 1 3906E+01 

.85278E+00 

-.20753E+00 

. 16969E-01 

ABSORPTION 

RMS 

Error  =0.6% 

\ 

0 

1 

2 

3 

0 

.63909E-01 

.48691 E-01 

-.11011 E-01 

.36647E-03 

1 

- . 25430E-01 

-.111 09E+00 

. 26393E-01 

-.211 54E-02 

2 

•11770E+00 

. 36211 E+00 

-.85429E-01 

.68359E-02 

3 

-.13761E+00 

-.181 36E+00 

•42026E-01 

-.33420E-02 

Table  10.  Approximating  predictor  coefficients,  b^. , 
for  a cirrus  cloud 


REFLECTION 

RMS  Error  =1.5% 

X- 

0 

1 

2 3 

0 

1 

2 


.25264E+00 
••29273E+0Q 
•. 24531 E+00 


. 45462E+00 
.22029E+01 
.63610E+01 


•29492E+00 
■24834E+01 
. 75864E+01 


.82536E-01 

.82674E+00 

•.25896E+01 


3 ! 

. 35097E+00 

. 38933E+01 

-.49188E+01 

.17187E+01 

TRANSMISSION 

RMS 

Error  =1.7% 

0 

1 

2 

3 

0 

.78722E+0Q 

-.66348E+00 

.44495E+00 

-. 12500E+00 

1 

2 


. 73255E+00 
. 68508E+00 


••43772E+01 

.11764E+02 


.47094E+01 
. 13716E+02 


■. 15549E+01 
.46796E+01 


3 

. 1 5453E+00 

-.70865E+01 

.8/472E+01 

-.30566E+01 

ABSORPTION 

RMS 

Error  = 0.1% 

0 

1 

2 

3 

0 

.49388E-01 

. 1 371 5E+00 

10122E+00 

.28402E-01 

1 

- . 1 9622E+00 

.11882E+01 

-.12267E+01 

. 40299E+00 

2 

. 241 1 8E+00 

-.23149E+01 

•28776E+01 

-.10168E+01 

3 

-.94420E-01 

. 1 1 305E+01 

-.15772E+01 

.58577E+00 
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PARAMETERIZATION  OF  THE  SURFACE  REFLECTION 

The  parameterization  of  reflection,  transmission  and  abosrption 
of  solar  radiation  by  clouds  described  in  the  previous  chapter  does 
not  include  the  effect  of  multiple  reflections  between  the  earth's 
albedo  surface  and  the  cloud  base.  Schneider  and  Dickinson  (1976)  em- 
phasized the  importance  of  these  multiple  reflections  in  determining 
the  solar  flux  reaching  the  surface,  especially  in  climate  models 
which  predict  snow  and  ice  cover. 

To  include  the  surface  reflection  and  the  atmospheric  effect  be- 
tween the  cloud  and  the  earth's  surface,  the  surface  is  assumed  to 
reflect  according  to  Lambert's  law  with  a surface  albedo  of  A$. 
Rayleigh  scattering  between  the  cloud  and  surface  is  not  considered 
and  the  absorption  between  the  cloud  and  surface  is  assumed  to  be 
mainly  due  to  water  vapor.  The  average  reflection,  transmission  and 
absorption  of  solar  radiation  by  cloud  layers  is  then  defined  as 
follows: 


1 


r - 2 / 

0 

r(p0>y0  dP0 

(17) 

1 

t - 2 / 

0 

d^o 

(18) 

1 

a » 2 / 

0 

a^oK  d^0 

(19) 
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Note  that  r is  equivalent  to  the  spherical  albedo  which  represents  the 
reflecting  power  of  the  entire  planet. 

Figure  13  depicts  the  surface  reflection  contribution  to  the  up- 
ward flux  at  the  cloud  base.  This  additional  fraction  of  the  upward 
flux  reaching  the  cloud  bottom  may  be  expressed  in  terms  of  an  infinite 
series  as  follows: 


f(y0) 


= t(u0)[(l-a)2As+(l-a)4As2  F + (1-a)6  A$3  r 2 + ...J 
- t(uQ)(l-a)2  As[l+(l-a)2  A$  r + (1-a)4  A$2  F2  + "*] 


(l-a)*  A 

- t(y  ) Hz 

0 l-(l-a)2  Asr 


(20) 


where  a represents  the  water  vapor  absorption  between  the  cloud  base 
and  the  earth's  surface.  Letting  the  superscript  * denote  parameters 
which  include  the  surface  effect,  the  actual  cloud  reflection,  trans- 
mission and  absorption  values  are  then  given  by 

r*(uQ)  = r(yQ)  + f(yQ)  t , (21) 

t*(uQ)  = t(yQ)  + f(yQ)  r , (22) 

a*(uQ)  * a(yQ)  + f(yQ)  a . (23) 

Equations  (21)  and  (22)  are  equivalent  to  those  derived  by  Chand- 
rasekhar (1950)  for  the  reflected  and  transmitted  intensities  In  the 
planetary  problem  based  on  the  principles  of  Invariance. 

The  water  vapor  absorption  a depends  on  the  water  vapor  path 
length  and  atmospheric  pressure.  Liou  and  Sasamori  (1975)  modified 


the  spectral  absorptivities  of  water  vapor  and  carbon  dioxide  for 
strong  and  weak  bands  based  on  measurements  of  Howard  et  al.  (1956) 
and  gave  a general  expression  in  the  form 
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where 


X1  = uP 


Ki/°i 


Xoi  = 10 


-Ci/Di 


(24) 


In  the  above  equations,  Av  is  the  spectral  width  of  the  absorption 
band,  u the  water  vapor  path  length,  P the  atmospheric  pressure  and 
C,  D and  K are  constants.  Thus,  the  total  water  vapor  absorption  be- 
tween the  cloud  and  the  surface  may  be  obtained  by 

a(u,F)  = l , (25) 

where  the  summation  is  for  all  the  water  vapor  absorption  bands  in  the 
solar  spectrum  and  P denotes  the  mean  pressure  between  the  cloud  and 
surface. 

Equations  (21)  - (23)  were  used  to  compute  the  reflection,  trans- 
mission and  absorption  values  for  non-zero  surface  albedos.  These 
values  were  compared  with  those  calculated  from  the  Inhomogeneous 
transfer  program  Involving  the  surface  reflection.  Figure  14  shows 
this  comparison  for  a cumulus  cloud  In  a tropical  atmosphere  with  a 
surface  albedo  of  0.4.  Comparisons  were  made  for  different  cloud  types, 
atmospheres,  and  surface  albedos  and  revealed  that  differences  between 
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the  two  methods  were  generally  less  than  about  2%. 

Calculations  were  also  made  using  the  surface  reflection  parame- 
terization as  described  previously  except  that  the  water  vapor  absorp- 
tion between  the  cloud  base  and  the  earth's  surface  was  neglected; 
i.e.  a is  zero  in  Eq.  (20).  The  values  for  reflection  and  transmission 
showed  no  significant  difference  from  those  which  considered  the  ab- 
sorption by  water  vapor;  however,  the  error  in  absorption  values  was 
about  2-3%  greater. 

It  was  shown  in  the  previous  chapter  that  the  error  of  the  parame- 
terization of  the  solar  radiative  properties  of  clouds  was  less  than 
2%  for  each  of- the  four  cloud  types.  In  addition.  Figure  14  shows  that 
the  error  of  the  surface  reflection  parameterization  is  also  less  than 
2%.  Figure  15  displays  the  comparison  between  reflection,  transmis- 
sion and  absorption  values  calculated  from  the  inhomogeneous  transfer 
program  and  from  the  combination  of  the  two  parameterizations  described 
previously.  These  values  are  again  for  a cumulus  cloud  in  a tropical 
atmosphere  with  a surface  albedo  of  0.4.  Comparisons  made  with  the 
four  cloud  types  showed  that  the  difference  of  values  was  normally  less 
than  3%. 
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CHAPTER  5 
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PARAMETERIZATION  OF  REFLECTIVITY,  TRANSMISSIVITY  AND 
EMISSIVITY  OF  INFRARED  RADIATION  BY  A CIRRUS  CLOUD 

Since  cirrus  clouds  may  not  be  considered  black  bodies  as  are 
other  cloud  types,  their  infrared  radiative  properties  are  quite  impor- 
tant in  studies  of  the  weather,  climate  and  radiation  balance  of  the 
earth-atmosphere  system.  The  infrared  portion  of  the  radiative  trans- 
fer model  was  used  to  obtain  the  fluxes  necessary  to  compute  the  reflec- 
tivity, transmissivity  and  emissivity  of  infrared  radiation  by  a cirrus 
cloud.  These  infrared  radiative  properties  Of  a cirrus  cloud  were  com- 
puted over  the  range  of  vertical  ice  contents  listed  in  Table  6.  The 
flux  of  radiation  emitted  by  the  atmosphere  and  entering  the  cloud  base, 
F (t^)  and  the  flux  of  radiation  due  to  atmospheric  emission  entering 
the  cloud  top,  F'*'^) , were  set  equal  to  zero  when  calculating  the 
emissivity  of  the  cloud  so  that  the  effect  of  atmospheric  emission  was 
not  included.  On  the  other  hand,  the  Planck  function  within  the  cloud 
layer  was  set  equal  to  zero  when  calculating  the  reflectivity  and 
transmissivity  In  order  to  exclude  the  effect  of  the  cloud  emission. 

The  plots  of  the  computed  radiative  properties  are  displayed  in 

Figure  16.  From  these  plots  it  can  be  seen  that  each  of  the  radiative 

properties  become  nearly  constant  beyond  a vertical  ice  content  of 

100  g • m (approximately  2 km  geometrical  thickness);  therefore,  com- 

-2 

putations  beyond  150  g • m were  not  required.  There  is  also  a 
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Figure  16.  Infrared  radiative  properties  of  a cirrus  cloud. 
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significant  difference  between  the  emissivities  of  the  cloud  at  its 
base  and  top.  This  is  due  to  the  temperature  gradient  within  the  cloud 
layer  (T^  > T ).  Since  the  emission  of  the  cloud  is  dependent  upon  its 
temperature,  the  emission  of  the  cloud  sub-layers  below  the  cloud  top 
will  be  greater  than  the  emission  of  the  sub-layer  at  the  cloud  top. 
This  will  cause  a positive  flux  divergence  in  the  upward  direction  and 
a negative  flux  divergence  in  the  downward  direction.  Consequently, 
the  cloud  emissivity  is  less  at  the  cloud  base  than  at  the  cloud  top. 

The  reflectivity,  transmissivity  and  emissivity  curves  were  fit 
to  an  approximating  polynomial  equation  using  a least  squares  technique 
similar  to  that  described  in  Chapter  3.  The  following  polynomial  equa- 
tion provided  the  best  fit  of  the  data: 


R(W)  = 


(26) 


where  R(W)  denotes  the  flux  reflecti vity,  transmissivity  or  emissivity 

of  the  cirrus  cloud,  c^  the  predictor  coefficients  and  W the  vertical 
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ice  content  in  units  of  10  g-m  . 

Table  11  lists  the  coefficients  which  were  computed  for  the  vari- 
ous infrared  radiative  properties  of  the  cirrus  cloud.  In  each  case, 
the  root  mean  square  error  of  the  approximating  functions  is  less  than 
IX. 


Table  11.  Approximating  predictor  coefficients,  c^ , 
for  a cirrus  cloud 


i 

Reflectivity 

Transmissivity 

Emissivity 

(Cloud  Top)  (Cloud  Base) 

0 

. 30619E-01 

.73597E+00 

'.28363E+00 

.28773E+00 

1 

.81 1 34E+00 

-.38162E+01 

•36135E+01 

. 34309E+01 

2 

-.18995E+01 

.83288E+01 

-.77273E+01 

-.75277E+01 

3 

.24900E+01 

- .91066E+01 

.83695E+01 

.82498E+01 

4 

- . 1 5805E+01 

.48765E+01 

-.44646E+01 

-.44163E+01 

5 

. 37531 E+00 

-. 10141E+01 

.92724E+00 

.91693E+00 

CHAPTER  6 


THE  RADIATION  BALANCE  OF  THE  EARTH-ATMOSPHERE  SYSTEM 

The  parameterized  equations  which  were  derived  for  the  solar  and 
Infrared  radiative  properties  of  clouds  were  used  in  a relatively 
simple  solar  and  infrared  radiation  model  in  order  to  determine  the 
various  components  of  the  radiation  balance  of  the  earth-atmosphere 
system.  The  following  sections  of  this  chapter  include  descriptions  of 
the  solar  and  Infrared  portions  of  the  radiation  balance  model  and  the 
complete  radiation  balance  which  was  calculated. 

6.1  The  Global  Solar  Model 

The  solar  radiation  balance  of  the  earth-atmosphere  system  was 
calculated  using  the  transfer  program  for  the  cloudless  portions  of  the 
atmosphere  and  the  system  of  parameterized  equations  for  the  cloudy 
portions  of  the  atmosphere.  The  major  factors  considered  in  the  deter- 
mination of  the  radiation  balance  are  the  atmospheric  profile,  the  geo- 
metrical properties  of  clouds,  the  fractional  cloudiness  for  each  cloud 
type,  the  earth's  surface  albedo,  the  duration  of  sunlight,  and  the 
zenith  angle  of  the  sun. 

The  radiative  transfer  program  calculations  for  the  cloudless  por- 
tion of  the  earth-atmosphere  system  utilized  atmospheric  profiles  of 
pressure,  temperature,  density,  water  vapor  and  ozone  which  were 
derived  by  averaging  the  five  sets  of  atmospheric  profiles  described 
In  Chapter  2. 
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The  cloud  geometrical  properties  which  are  of  interest  are  height 
in  the  atmosphere,  thickness  and  horizontal  extent.  Table  12  lists 
the  mean  values  which  were  obtained  by  averaging  the  cloud  height  dis- 
tribution values  given  by  Telegadas  and  London  (1954)  and  the  distri- 
butions of  fractional  cloudiness  given  by  London  (1957)  for  the  North- 
ern Hemisphere  and  by  Sasamori  et  al.  (1972)  for  the  Southern  Hemi- 
sphere. These  average  values  take  into  account  the  variations  with 
season  and  latitude.  In  addition,  the  fractional  cloudiness  values 
for  nimbostratus  and  cumulonimbus  clouds  were  included  in  the  frac- 
tional cloudiness  of  the  altostratus.  All  three  cloud  types  have  rela- 
tively large  optical  depths  and  therefore  the  solar  radiative  proper- 
ties are  quite  similar. 


Table  12.  Mean  Cloud  Heights  and  Global  Fractional  Cloudiness 


Cloud  Type 

Cloud  Base 
(km) 

Cloud  Thickness 
(km) 

Fractional 

Cloudiness 

Cumulus 

1.55 

0.71 

0.09 

Altostratus 

2.96 

0.80 

0.18 

Stratus 

1.18 

0.10 

0.12 

Cirrus 

8.35 

1.70 

0.14 

The  surface  albedo  of  the  earth  Is  also  an  Important  parameter 
since  It  determines  the  amount  of  transmitted  solar  radiation  reaching 
the  surface  which  Is  reflected  back  into  the  atmosphere  to  be  absorbed 
or  scattered,  or  to  escape  back  into  space  as  a component  of  the 
earth's  global  albedo.  Surface  albedo  values  for  the  Northern 
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Hemisphere  were  taken  from  the  work  of  Katayama  (1967b),  and  for  the 
Southern  Hemisphere  from  Sasamori  et  al . (1972).  These  values  were 
averaged  with  respect  to  latitude  and  season.  The  global  surface 
albedo  was  calculated  to  be  15%. 

The  duration  of  sunlight  and  the  solar  zenith  angle  are  also  im- 
portant parameters  in  determining  the  radiation  balance  of  the  earth- 
atmosphere  system.  The  zenith  angle  of  the  sun  is  computed  from 

sina  = sin$  sin6  + cos<J>  cos6  cosh  (27) 

where  a is  the  altitude  of  the  sun  (angular  elevation  above  the  hori- 
zon), <p  the  latitude  of  the  observer,  5 the  declination  of  the  sun, 
and  h the  hour  angle  of  the  sun.  The  cosine  of  the  solar  zenith  angle 
is  then  given  by 

yQ  ■ cose  * cos(90°-a)  (28) 

In  general,  the  solar  zenith  angle  varies  significantly  each  hour  of 
the  day  and  sunlight  duration  varies  with  season  as  well  as  latitude. 
Table  13  shows  the  weighting  factors  which  were  calculated  for  the  six 
values  of  used  In  the  transfer  program.  These  weighting  factors 
compensate  for  the  variations  In  solar  zenith  angle  and  sunlight 
duration. 


Table  13.  Weighting  Factors  for  the  Cosine  of  Solar  Zenith  Angle 
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6.2  The  Global  Infrared  Model 

The  global  infrared  radiation  budget  for  the  earth-atmosphere 
system  is  somewhat  easier  to  model  than  the  solar  radiation  budget 
since  there  is  no  dependence  on  surface  albedo  or  zenith  angle.  There 
fore,  the  factors  which  must  be  considered  in  order  to  compensate  for 
the  seasonal  and  latitudinal  variations  are  the  atmospheric  profile, 
the  geometrical  properties  of  clouds,  and  the  fractional  cloudiness  of 
each  cloud  type. 

The  infrared  portion  of  the  global  radiation  balance  was  calcu- 
lated using  both  the  transfer  program  and  the  parameterized  equations 
.for  the  infrared  radiative  properties  of  a cirrus  cloud.  From  the 
transfer  program,  values  for  the  upward  fluxes  at  the  atmospheric  top 
and  the  downwelling  radiation  at  the  earth  surface  were  obtained  using 
the  averaged  atmospheric  profile,  described  In  the  preceding  section, - 
and  the  averaged  values  of  cloud  distribution  listed  In  Table  12.  The 
upward  flux  at  the  surface  was  computed  using  the  surface  temperature 
in  the  Stefan-Boltzmann  law 

F = aTs4  (29) 

In  this  equation,  F is  the  total  flux  of  infrared  radiation  emitted  by 
the  earth  and  a is  the  Stefan-Boltzmann  constant  which  has  a value  of 
8.128  x 10"11  cal  cm-2  min"1  V4. 

6.3  The  Global  Radiation  Balance 

The  annual  radiation  budget  of  the  earth-atmosphere  system  is 
presented  In  Figure  17.  The  radiation  from  the  sun  on  the  earth,  aver 
aged  for  the  entire  year,  is  represented  by  100  units.  Using  a solar 


Figure  17.  Radiation  balance  of  the  earth-atmosphere  system. 
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constant  of  1353  W-m  (Thekaekara,  1976),  the  average  insolation  at 

-2  -1 

the  top  of  the  atmosphere  was  calculated  to  be  0.485  cal -cm  -min  . 
Basically,  the  figure  is  in  three  sections,  one  dealing  with  solar 
radiation  and  the  manner  in  which  it  is  apportioned  in  the  atmosphere; 
the  second  concerning  infrared  radiation  and  its  distribution;  and  the 
third  dealing  with  nonradiative  processes.  For  the  solar  part,  the 
important  parameters  are  the  absorption  of  solar  radiation  by  the 
cloudless  atmosphere  and  by  clouds,  absorption  and  reflection  by  the 
earth's  atmosphere,  and  the  various  components  of  the  reflected  radia- 
tion at  the  top  of  the  atmosphere.  The  parameters  of  Interest  for  the 
infrared  radiation  include  the  upward  flux  at  the  earth's  surface,  the 
upward  flux  at  the  atmospheric  top,  and  the  net  loss  of  Infrared  radia- 
tion from  the  atmosphere.  The  net  upward  flux  at  the  surface  is  the 
difference  between  the  flux  emitted  by  the  surface  and  the  downward 
flux  from  the  atmosphere  reaching  the  surface.  This  Is  always  a posi- 
tive quantity  since  the  black-body  emission  from  the  earth  is  always 
greater  than  the  non-black  emission  from  the  cooler  atmosphere.  The 
net  loss  from  the  atmosphere  is  the  difference  between  the  upward  flux 
at  the  top  of  the  atmosphere  and  the  net  upward  flux  at  the  surface. 

The  radiative  balance  between  the  surface  of  the  earth  and  the 
atmosphere  Involves  fluxes  of  latent  and  sensible  heat  In  addition  to 
the  radlatlonal  fluxes.  In  this  study,  no  attempt  has  been  made  to 
make  an  Independent  estimate  of  these  nonradiative  components  of  the 
radiation  balance.  The  flux  of  latent  and  sensible  heat  represented 
In  Figure  17  has  been  selected  to  produce  an  overall  balance  at  the 
surface.  Consequently,  it  Includes  the  Integrated  errors  of  the  other 
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components  of  the  global  radiation  balance.  The  average  annual  ratio 
of  sensible  to  latent  heat  loss  at  the  surface  (Bowen  ratio)  is  approx- 
imately 0.22  for  the  Southern  Hemisphere  and  0.32  for  the  Northern 
Hemisphere  (Sasamori  et  al.,  1972).  Therefore,  the  global  value  of  the 
latent  heat  component  is  23  units  and  the  sensible  heat  component  is 
6 units. 

In  Table  14  comparisons  between  the  present  study  and  previous 
models  are  presented  for  the  solar  radiation  absorbed  and  reflected  by 
the  earth-atmosphere  system,  the  net  upwelling  infrared  radiation  at 
the  earth's  surface,  the  total  upwelling  infrared  radiation  at  the 
atmospheric  top,  the  net  infrared  radiation  absorbed  by  the  atmosphere, 
and  the  latent  and  sensible  heat  transport.  It  should  be  noted  that 
the  works  of  Houghton  (1954)  and  London  (1957)  are  for  the  Northern 
Hemisphere  and  that  of  Sasamori  et  al . (1972)  is  for  the  Southern  Hemi- 
sphere only. 

The  differences  in  the  absorption  and  reflection  of  solar  radia- 
tion by  clouds  between  the  various  studies  Is  partially  due  to  the 
different  cloud  distributions,  but  is  primarily  caused  by  the  parame- 
terization methods  used  in  the  models.  Houghton  used  a single  set  of 
cloud  albedos  for  the  entire  model,  while  London  and  Sasamori  used  a 
range  of  values  dependent  upon  cloud  thickness  and  to  some  extent  upon 
the  solar  zenith  angle,  but  do  not  Include  the  additional  effect  of 
the  surface  reflection  which  Is  Included  In  the  parameterization 
developed  In  this  study. 

The  total  absorption  of  solar  radiation  computed  In  the  present 
work  exceeds  the  absorption  computed  by  the  previous  investigators. 


Table  14.  Annual  Radiation  Budget  of  the  Earth-Atmosphere  System 
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This  is  due  to  the  increased  absorption  from  aerosols  and  the  effects 
of  scattering  by  aerosols  and  clouds,  which  increase  the  amount  of 
energy  available  for  absorption  by  increasing  the  optical  path  length 
through  the  atmosphere. 

The  comparisons  of  planetary  albedo  show  that  the  present  calcu- 
lations yield  an  albedo  which  is  3-4%  lower  than  those  of  the  older 
works;  however,  the  planetary  albedo  value,  0.31,  of  this  work  is  in 
good  agreement  with  those  derived  from  satellite  observations.  Vonder 
Haar  and  Suomi  (1971)  calculated  an  albedo  of  0.30  from  TIROS,  Nimbus 
and  ESSA  satellite  measurements  for  the  period  1962-1966.  Raschke  and 
Bandeen  (1970)  measured  an  albedo  of  0.29  to  0.31  for  June  and  July  of 
1966  from  Nimbus  2 satellite  observations.  Raschke  et  al.  (1973), 
using  Nimbus  3 measurements  from  1969-1970,  evaluated  the  global  albedo 
at  0.284.  More  recently.  Smith  et  al . (1977),  using  Nimbus  6 earth 
radiation  budget  experiment  data,  computed  the  global  albedo  to  be 
0.30.  The  difference  of  the  planetary  albedo  derived  by  this  model  and 
the  satellite  observed  values  are  probably  due  to  the  overestimation 
of  cloudiness  and  the  underestimation  of  absorption  in  the  atmosphere 
due  to  the  uncertainties  and  Inaccuracies  in  the  treatment  of  aerosols. 

The  Infrared  radiation  emitted  from  the  earth's  surface  and  the 
atmosphere,  as  computed  by  the  model  used  In  this  study,  are  somewhat 
greater  than  those  calculated  previously.  Consequently,  the  net  radi- 
ation from  the  earth's  surface  is  smaller  and  the  net  radiation  lost 
by  the  atmosphere  Is  greater  than  those  of  the  earlier  works.  These 
differences  result  mainly  from  the  different  water  vapor,  ozone,  tem- 
perature and  cloud  distributions  used.  The  parameterization  of  the 
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radiative  properties  of  the  non-black  cirrus  cloud  also  plays  a signif- 
i icant  role  in  the  overall  emission  of  the  earth-atmosphere  system. 

. The  present  study  is,  however,  in  good  agreement  with  the  satellite 

1 observed  upwelling  infrared  flux.  The  studies  of  Smith  et  al . (1977) 

and  Vonder  Haar  and  Suomi  (1971)  calculated  values  which  equal  70%  of 
the  incoming  solar  radiation. 
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CHAPTER  7 

CONCLUSIONS 

Clouds  are  the  most  important  atmospheric  elements  involved  in 
the  complex  interactions  of  the  radiation  field.  In  this  study,  the 
cloud  effects  were  carefully  reproduced  by  the  comprehensive  radiative 
transfer  program  and  then  parameterized  into  simplifying  empirical- 
theoretical  equations  for  both  the  solar  and  infrared  radiative  proper- 
ties of  clouds.  The  reflection,  transmission  and  absorption  of  solar 
radiation  by  four  cloud  types  are  parameterized  by  polynomial  equa- 
tions in  terms  of  the  solar  zenith  angle  and  the  cloud  vertical  liquid 
water/ice  content.  Parameterized  equations  for  the  infrared  reflec- 
tivity, transmissivity  and  emissivity  of  cirrus  clouds  are  expressed 
in  terms  of  the  cloud  vertical  ice  content. 

By  employing  various  different  atmospheric  profiles,  the  effect  of 
the  atmospheric  profile  on  these  radiative  properties  of  clouds  is 
shown  to  be  insignificant.  However,  the  cloud  type,  which  is  repre- 
sented by  a distinct  particle  size  distribution,  does  have  a signifi- 
cant Influence  on  the  values.  The  additional  effect  of  the  surface 
reflection  upon  the  reflection,  transmission  and  absorption  of  solar 
radiation  by  clouds  is  further  parameterized  In  which  water  vapor 
absorption  between  the  cloud  and  the  earth's  surface  is  taken  into 
account. 

The  comprehensive  radiative  transfer  program  is  also  utilized  in  a 
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a relatively  simple  radiation  budget  model.  The  global  results  ob- 
tained are  in  good  agreement  with  earlier  radiation  budget  studies  and 
agree  more  closely  with  recent  satellite  observations  than  do  the 
earlier  works. 

Overall,  it  is  felt  that  the  parameterized  equations  for  the 
solar  and  infrared  radiative  properties  of  clouds  proved  accurate  in 
the  radiation  budget  model  and  would  be  useful  and  practicable  in  con- 
nection with  the  study  of  the  climate  and  climatic  changes  of  the 
earth-atmosphere  system. 
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